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Abstract 

We discuss recent measurements of the wavelength-dependent absorption coefficients in 
deep South Pole ice. The method uses transit time distributions of pulses from a variable- 
frequency laser sent between emitters and receivers embedded in the ice. At depths of 800 
to 1000 m scattering is dominated by residual air bubbles, whereas absorption occurs both 
in ice itself and in insoluble impurities. The absorption coefficient increases approximately 
exponentially with wavelength in the measured interval 410 to 610 nm. At the shortest 
wavelength our value is about a factor 20 below previous values obtained for laboratory 
ice and lake ice; with increasing wavelength the discrepancy with previous measurements 
decreases. At ~ 415 to ~ 500 nm the experimental uncertainties are small enough for us to 
resolve an extrinsic contribution to absorption in ice: submicron dust particles contribute by 
an amount that increases with depth and corresponds well with the expected increase seen 
near the Last Glacial Maximum in Vostok and Dome C ice cores. The laser pulse method 
allows remote mapping of gross structure in dust concentration as a function of depth in 
glacial ice. 
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1 Introduction 



The AMANDA experiment (Antarctic Muon and Neutrino Detector Array) is a 
U.S. -Germany-Sweden collaboration with the aim of deploying a facility to detect 
neutrinos from astrophysical sources (a "neutrino telescope" ) deep in the South Pole 
ice cap Q. High-energy muon neutrinos will convert into charged muons through 
charged weak current interactions somewhere in the ice volume (or in the bedrock 
below the ice cap). As the muons are very penetrating, and move at relativistic 
speed, they will generate Cherenkov light along, typically, several hundred meters 
of track length. By recording the arrival times of Cherenkov photons at several 
locations in the ice volume, the direction of the muon track can be reconstructed. 
Interesting events will be muons coming from below the horizon, since they can only 
have been generated by upcoming neutrinos (all other known particles would be 
stopped in the interior of the earth). 

The AMANDA detector will eventually consist of many strings, each with 20-50 
optical modules containing photomultiplier tubes (PMTs) lowered around 2 kilome- 
ters into the ice. Typical spacing between strings and between modules in strings 
will be tens of meters. 

The performance of the detector will depend crucially on the optical properties 
of the ice. A first stage of the detector, deployed during the 1993-94 austral summer, 
consists of four vertical strings each having 20 PMTs spaced at 10 m distance. The 
strings form an 80 deg sector of a circle of 30.4 m radius, with one string at the center 
and three strings on the circumference. A laser calibration system was developed 
in order to determine the precise detector geometry and the optical properties of 
the ice. A brief summary of the results has recently been published in [§, 9. The 
results have shown that the ice at 800 - 1000 m depth is highly turbid. The short 
scattering lengths (between 10 and 20 cm) are consistent with being due to a residual 
population of bubbles of trapped air. Deeper in the ice these bubbles are expected 
to be transformed to air hydrate crystals, which have a refractive index close to that 
of ice § and therefore should not cause much scattering. (For a thorough discussion 
of various scattering mechanisms in deep Antarctic ice, see the accompanying paper 

Surprisingly, the absorption coefficient of the ice at wavelengths between 400 and 
600 nm was found to be smaller by a large factor (greater than 10 at the shortest 
wavelengths) than the value for laboratory ice pT| , ^] (for a compilation of several 
sets of experimental data see ||). The South Pole ice is thus a remarkably clear 
optical medium, and in the 1995-96 season we have deployed 4 new strings to a 
maximum depth of 2000 m, where air bubbles should have entirely transformed into 
air hydrates. The analysis of data from this deployment is in progress. 

Since the reported AMANDA results (especially those of ||) deviate from pre- 
viously published values by such a large factor, we have felt the need to present 
to the optical physics community the details of our experimental method, the pro- 
posed explanation and implications of the long absorption lengths and an analysis 
of absorption in dust present in the ice. Various sources of residual scattering at 
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bubble- free depths are discussed in ||. 

2 Outline of the method to extract scattering and ab- 
sorption lengths 

The experimental setup for the measurement of the optical properties of the South 
Pole deep ice is shown in Fig. 1. Light from a dye laser (driven by a Nd:YAG 
laser) was brought down optical fibers to diffusing nylon spheres (one per PMT) 
approximately 30 cm below each PMT. 

Pulsed laser light (4 ns FWHM) was fed into the optical fibers at a frequency 
of 10 Hz. By using gratings in the dye laser and a monochromator the wavelength 
of the light could be varied between 410 and 610 nm, with a precision of around 2 
nm. (At shorter wavelengths, not enough light could be transmitted through the 
kilometer-long fibers.) After passage down the fiber, the 4 ns pulses had broadened 
to around 12 ns. Light was led down to one nylon sphere at a time, and the response 
of neighboring PMTs in the same string and in the other strings was recorded. 

The PMTs (8 inch EMI 9353 tubes with a noise rate in the ice of about 1.7 kHz) 
were operated at a gain of 10 8 , and had approximately 6 ns resolution, including 
time smearing in the coaxial cables that fed the signals up to the surface. The signal 
from the PMT 30 cm from the flashing nylon sphere defined t=0 (corrected for run 
time through fibers, etc.). 

The main purpose of the laser system just described was to serve as a calibration 
system for the absolute timings of the PMTs and their relative positions in the 
AMANDA muon detector. Immediately after the system became operational it 
was realized, however, that the South Pole ice at depths between 800 and 1000 
meters is not an ideal optical medium for imaging the tracks of muons through their 
Cherenkov light. Instead of a narrow arrival pulse (expected, typically, at 87 ns for 
20 m distance between emitting and receiving tube), a very broad structure with 
maximum at around 500 ns was observed (see Fig. 2). As very exact timings are 
required to recontruct the direction of muons from the sparsely sampled Cherenkov 
photons they generate, the observed timing pattern makes the present array not 
useful for muon tracking (but the array is used for other astrophysical applications 

However, the presence of the broad time distributions, and our modeling of them, 
have permitted us to make more accurate measurements of the optical properties 
of ice in the visible region than have ever been done in the laboratory. As shown 
below, the in situ laser measurements in large volumes of ice and the random walk 
nature of the light propagation permit a clear separation between absorption and 
scattering, something which is extremely difficult to study in the laboratory. 

In that process, we have independently discovered a method similar to one re- 
cently used in medical and other applications pi], one difference being the longer 
time scales and distances involved here, another the fact that both the light emit- 
ters (nylon spheres) and receivers (PMT tubes) are immersed in the medium to be 
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measured. The derivation that we will present of the relevant formulas seems to be 
new (and simple), and may be of interest for optical physics researchers. We will 
also show that a detailed Monte Carlo simulation of light propagation supports the 
analytical treatment, and fits the measured distributions excellently. 



2.1 Random walk description of light propagation 

The physical picture of photon propagation in the bubbly ice is that of a random 
walk. When a photon leaves the nylon sphere it passes through clear ice along 
a straight trajectory for some path length until it hits an air bubble. Then it is 
refracted (or externally reflected), travelling in a new direction until it hits the next 
bubble. The wide dispersion of arrival times at the PMTs must mean that a typical 
photon has undergone many scatterings, and therefore a probabilistic analysis of the 
random walk should be possible. The problem is simplified by the randomness of 
the medium and the macroscopic distances of propagation which means that any 
interference effects can be ignored. 

We define the simple, geometrical scattering length on spherical bubbles as 

A(m6 = , (1) 

C-bub 

with the scattering coefficient c bub given by 

f°° dn bub (r) 2 

Cbub = / 3 Trr dr, (2) 

Jo dr 

where dn bub (r)/dr is the number density of bubbles having radius between r and r + 
dr; itr 2 is the geometrical cross section. This means that the probability distribution 
for the distance d between two successive scatterings is proportional to e~ d ^ bub . We 
note that the average geometrical scattering length is 

(d) = \ bub (3) 

and the second moment is 

(d 2 ) = 2AL (4) 

We only have to incorporate two modifications to the standard random walk 
treatments given in the literature (e.g, [||), namely absorption of the light beam 
and the non-isotropy of the scattering amplitude for photons incident on air bub- 
bles. Since the air bubbles themselves do not cause any absorption of photons, only 
scattering, it is very simple to incorporate absorption in the random walk picture. 
Suppose that a photon arrives at a PMT after N steps. Then the total path length 
is on the average NX bub , and absorption can be taken into account by weighting that 
particular path by the factor 

~ N *bub , , 

e ^ a (5) 
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(of course we assume N >> 1; typically ./V is on the order of 1000 in our problem), 
where A a is the true absorption length (i.e. related to true absorption of photons in 
the medium). 

Sometimes it is convenient to discuss absorption in terms of the absorption co- 
efficient, 

O-abs = t — • (6) 

■^abs 

This is related to the absorptive (imaginary) part of the index of refraction by 

47rlm(n) 

Oiabs = 7 , (') 

where A^ is the wavelength of the light source. 

Non-isotropy of the scattering from bubbles means that there is a correlation 
between successive vectors making up the random walk. If the bubbles are spherical 
there is still an azimuthal symmetry, and one can show that (e.g. by introducing 
rotation matrices which rotate each successive vector to the polar axis) 

in-n+k) = aLOW + A (8) 

where r = (cos 9) is the average of the cosine of the scattering angle, and the factor 
of 2 entering for k = (<5o,fc is a Kronecker delta function) is a reflection of Eq.(^) 
and Eq.d§. 

In the formula for random walk with absorption introduced as just explained: 

1 



W N (R) = - -e ^ 2 )Ne~^~, (9) 

(27^)^/3)3/2 

it just remains to evaluate: 

N 

(R 2 )n = {(^) 2 )- (10) 

1=1 

By expanding the squared sum in ((X)i=i ?i) 2 ) one finds N terms with relative 
distance (namely, |f*i| 2 + |r*2| 2 + . . . 4- |r/v| 2 ), 2(N — 1) terms of relative distance 1 
(r*i • 7*2 + 7*2 • T\ + . . . +r)v • r/v-i), and in general 2(iV — k) terms of relative distance 
k. 

Using then Eq. (|8j) , we find that the factor 2A 2 ub can be taken outside and we 
arrive at 



N-l 



(R 2 ) N = 2X 2 bub £ (AT - k)r k 



11) 



k=0 



2iVAL 



(1 - r) + (t n+1 - t)/N 



1-t) 



In the limit of large N we obtain 
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(i?)iv = 2^,_ (12) 

With the identification NAt, u b = c{t, where q = c/n is the velocity of light in the 
ice, we obtain the formula corresponding to the Green's function for the radiative 
transport of a spherically symmetric laser pulse emitted at a distance d = at time 
t = 

^ ) = (4^3/2 ^^ ( 13 ) 

with d 2 = R 2 , and where the constant of diffusion D is given by 

D = ^£ (14) 
with the effective scattering length A e // related to through the formula 

X eff = ^L. (15) 

We note that in the expression for D only the refractive scattering part of the op- 
tical amplitude needs to be taken into account. If diffraction or any other small-angle 
scattering effect is included, making the geometrical bubble cross section increase 
from a to {B + l)a, then 

hub -» W(l + B), r -» (B + t)/(S + 1), 

and -D is easily seen to be invariant under these transformations. This has the 
simple physical interpretation that the diffusion power of the very forward-peaked 
diffraction has to be negligible compared to the large-angle scattering induced by 
refraction on bubble walls. 

The formula (|l^) can be put in a more convenient form as follows using two 
dimensionless quantities y and k. Defining y = t/y/J with 

/ = (16) 



and the dimensionless quantity 



one gets 



V3d 
2 VhXff 

C 



(17) 



u{f,y) = ^e~<y +l M (18) 

y 2 

with 

C= r . (19) 

/ 1 /4(4vrc J Ae///3)5 
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Using the basic integrals 



dyy- l / 2 e-<y +1 M = J dyy-^e^^ = ^ V 2k (20) 
one finds the average value of arrival time t for a given distance d 



i V A eff 



i.e. a linear relation. 

In fact, defining the general class of integrals 



fl(K,n)=l dyy- 3 * +n e-<y +1 M ( 22 ) 



one finds a very convenient recursion relation by taking the derivative of Eq.(|2_ 
with respect to k: 

Mk, n + 1) = ~fi(n, n) - //(«, n - 1), (23) 
with initial conditions for the recursion given by 

/j(«,0) = / z («,l) = \f^e~ 2K (24) 

(the equality of fj(n,0) and //(«, 1) , or more generally, //(«;, n) and //(«;, 1 — n) is 
easily demonstrated by letting y — > 1/y in the integral.) 

Using the arrival time distribution expressed in this scaling form, one can easily 
calculate all moments of the distribution, e.g., the expectation value 

E(y) = 1, (25) 

and the variance 

V( y ) = E(y*)-E(yf = ±. (26) 

2k 

The scaling form of the distribution was used in the numerical fitting procedure for 
the arrival times of the laser pulses. 

The value of time t pea k at which the distribution peaks is given by 



^--t^yh^-^-)- (27> 

For the case at hand, with A a >> A e //, this also gives an essentially linear depen- 
dence of ipeafc on d (compared to t pea k oc d 2 for a random walk process without 
attenuation) . 
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2.2 Application to muons traversing bubbly ice 

The results for the Green's function derived here can also be used to obtain an 
expression for the time distribution of Cherenkov photons coming from a straight 
muon track, for which the closest approach ("impact parameter") di to the PMT at 
t = is much larger than the scattering length \ e ff (i- e - i n our case greater than 
a meter or so). Introducing a = 4X e ffCi/3 and a cut-off length I (the "size" of the 
detector) one finds 



f\Zc 3 /ci(t+l/c ) dz , , 

^ 6 I " --( v /cIci(z+l/z)+C2) 



I(di,t;l) = NJ^ -wcic 3 ^ /z)+ c 2) ^ (2g) 

V Ci Jmax^y/ca/ait-l/co)) Z A l l 

where N is a normalization factor, cq is the velocity of the muon (close to the velocity 
of light in vacuum) , and 

ci = (d? + c 2 t 2 )/a, (29) 

c 2 = -2c 2 t/a, (30) 

C3 = cl/a + Ci/\ a . (31) 

The total flux F^{di) received at a PMT of area Apmt from a muon generating 
Nc photons per unit track length passing the PMT at impact parameter d can be 
obtained (for di » Xbub) by integrating this expression over time and taking the 
limit I — ► oo. We find 



where i^o is a modified Bessel function of the second kind. 

In applications where directional information and time structure at the sub- 
microsecond level are not important, such as the detection of neutrinos from super- 
novae where relevant time scales are on the order of seconds, the present AMANDA 
detector in bubbly ice is perfectly functional. A neutrino-induced event from a SN 
explosion typically generates a positron track of only 10 cm, and can on the diffusion 
scale y/\ a X e ff be regarded as pointlike. Then the total response function (flux per 
unit area integrated over time) for N$ n photons (typically of the order of 3000 per 
positron) injected at distance d is given by 



isN(d) = j^-^y^. (33) 

We can define an effective volume of ice seen by a PMT by integrating Apmt 
times this flux over all space (the error caused by using the formula at small d where 
it is not strictly valid can be shown to be small, as it contributes a very small fraction 
of the result). 

We then obtain 
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t t N SN ApMT^a , 0/n 

vsn = t • [ay 

This result could (apart from the numerical factor) be anticipated on dimen- 
sional grounds using the observation that scatterings do not change the average 
photon flux (thus the result should not depend on the scattering length). For a 
real detector, correction factors related to the quantum efficiency of the PMTs and 
Poisson statistics of detected photons have to be added. Since the Cherenkov flux 
depends strongly on wavelength (~ 1/A 2 ) as do A a , the quantum efficiency and to 
some extent X e ff, an integration over wavelength has to be performed in addition. 

2.3 Monte Carlo simulation of light propagation in bubbly ice 

To verify that the general picture and derivations above are correct, we have also 
performed extensive Monte Carlo simulations of light propagation through the ice. 
Thanks to the existence of widely separated length scales, d ~ A a >> A e // » Xl, 
where A^ is the wavelength of the laser light, a fairly simple optical treatment is 
adequate to obtain an estimate of A a and X e ff good to a few percent accuracy or 
better. In fact, the estimate of A a is given quite directly by the fall-off of the arrival 
times of photons for large times, and is very insensitive to the uncertainties in the 
determination of A^. To determine X^ one needs the average value of the cosine 
of the scattering angle as input. One cause of uncertainty is that the bubbles may 
not be smooth enough to be regarded as perfect spheres. In that case, the diffusing 
ability could be greater (i.e. (cos 6) smaller) and A^;, consequently larger. Thus, we 
are able to determine X e ff quite accurately from laser data, but the relation between 
X e ff and the geometrical Xb u b contains uncertainties. A direct measurement of the 
scattering function of air bubbles contained in the ice cores taken from the Antarctica 
would obviously be of importance for this problem. 

Since photons from the coherent laser source scatter on the average between 
500 and 1500 times on bubbles before arriving at the PM tubes, the initial state 
in any single scattering can be regarded as being completely incoherent in phase, 
direction and polarization. We have thus employed simple geometrical optics for 
the scattering process incorporating singly refracted, externally reflected and totally 
reflected light with the relative strength of the components given by the Fresnel 
coefficients. (In fact, we have also performed full Mie scattering calculations and 
verified that the value of r obtained from this simple analysis is correct.) We take 
into account the change of index of refraction of the air in the bubble due to the 
pressure and temperature. For the ice we use the wavelength-dependent values of the 
refraction index given in Q. The bubbles are assumed to be spherically symmetric, 
in agreement with observations at these depths. 

Although a more complete treatment of the scattering could be implemented in 
the Monte Carlo, we emphasize that this, according to Eq. (|l5|), only amounts to 
a renormalization of the relation between the effective scattering length X e ff and 
the geometrical bubble scattering length Xb u b (with the latter defined as an average 
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over sizes and shapes of air inclusions), which is not important for discussion in this 
paper. (It is of importance, though, when discussing the implications of our results, 
e.g., for glaciology.) With our approximations we find for the quantity r entering 
Eq. (|i~5|), a value around 0.75. This is slightly wavelength-dependent as the refractive 
index of ice increases by some 2 % when going from 600 to 300 nm. The Cherenkov 
angle arccos(l/(nj/3)) (where (5 is the muon velocity in units of the speed of light 
in vacuum) will also increase by some 4 % in going from the longest to the shortest 
wavelengths. (Since the intensity spectrum of the Cherenkov radiation varies as 
1/A 2 , the shorter wavelengths are most important for muon detection. Below 300 
nm, however, the absorption in ice increases rapidly and the detection efficiency of 
the optical modules drops.) 



2.4 Fits to the experimental distributions 

In making global fits to the laser calibration data we used special routines that 
computed maximum likelihood fits to the observed distributions. The chi-squared 
method was used to give starting values for the maximum likelihood fit and to give 



a figure of merit of the fit. The Davidon variance metric algorithm [1C] was used for 
the minimization and the statistical errors were extracted from the variance matrix 
of the maximum likelihood procedure. The function used was formulated as 

,,\ x n • E n • Ax e /oc\ 

9{t) = 3— , (35) 

y^N 

where x n are the normalization constants for each measured arrival time distribution 
{x n ~ 1.0), E n the number of events in distribution n, Ax the bin size used in all 
distributions, 

T = Jy + - 

V y 
t t - At, 



y = —Ft or 



Vf Vf ' 

where At m allows for a shift in the definition of time zero for each emitting tube (this 
this depends on details of the PMT such as the value of the high voltage etc). We 
allowed for uncertainties in the vertical positions of the strings as that was difficult 
to determine during deployment. The initial chi-squared sum we minimized was 

X 2 = E E iVi ' n 9{ )) , (36) 

n =i i=i y^ n 

where the ith. bin in the nth distribution has a content of yi^ n events and the theo- 
retical prediction for this bin is g(t). 

There is a potential problem with this laser calibration method due to the fact 
that it is in general impossible to resolve several individual photons from the same 
laser pulse hitting the same PMT. Since the arrival time recorded by the receiving 
PMT is set by the first photon, one may tend to suppress long arrival times, which 
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could underestimate the absorption length. To overcome this problem, we moni- 
tored the total number of emitted laser pulses N. For each receiving tube we then 
estimated the average number of photons reaching this detector from the zero class 
of a Poisson distribution. If the distribution has n entries the Poisson parameter ji 
is given by 

P(0) = e-» = (37) 

"=K^y < 38 > 

The variable used to select non-saturated data is the sum of the zero class and the 
class with one photon, 

a = P(0) + P(l) = e' 1 * + ne'V = e" M (l + fj,), (39) 

i.e. 1 — q is the fraction of events expected to have two or more photons reaching the 
detector and thus causing saturation. To avoid the saturation problem, we demanded 
that 1 — a was less than 1 %. To get a statistically significant estimate of the optical 
parameters only time distributions where at least 2000 hits were registered by the 
receiving PM were selected. This left almost 500 time distributions for various 
wavelength and distance combinations. 

In Fig. 2 we show the distribution of arrival times as measured and as fitted to 
the analytical expression (Eq. 13) for our random walk model, for an emitter and a 
receiver separated by 10 m. The results for wavelengths 410 and 610 nm are shown 
as open and filled circles, respectively. The fits shown give the same value for the 
scattering length Xb u b (around 10 cm) for the two cases, but absorption lengths of 
around 230 m and 10 m, respectively. The quality of the analytical fit to the data 
is quite remarkable. We thus conclude that the physics behind the laser calibration 
data is well understood, and that we have, in fact, at our disposal a very powerful 
tool to make in situ measurements of important properties of the ice such as A a and 
Xbub- In particular, A a and A e // are essentially uniquely determined by the temporal 
distribution of arrived photons at each particular PM tube, once the geometry of 
the detector array is fixed. 

One may notice that the values of A a and X e ff obtained through this analy- 
sis are essentially spatial averages over the diffusion scale y/X a X e ff. Since we let 
these parameters depend on depth in the analysis we are using an "adiabatic ap- 
proximation" that should be excellent if the quantities vary slowly with depth. The 
goodness of the fit to the function ( |i"3| ) for all wavelengths and all depths shows the 
self-consistency of this scheme. 

3 Summary of Experimental Results 

In this Section, we briefly summarize the results on absorption length as a function of 
wavelength obtained by the AMANDA calibration system at two depths, around 830 
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and 970 m depth, respectively. In |j] results are also given for an intermediate region 
around 910 m, with results similar to those at the shallowest depth. The results are 
shown in Table 1 (for details, see ||]). The errors given are only statistical; a 
systematic error estimated to be around 5 % (dominated mainly by the uncertainty 
of the exact geometry of the AMANDA array should be added. 



4 Modeling the Optical Properties of Highly Pure South 
Pole Ice 

Using the results presented in Section 3, we next analyze the intrinsic optical prop- 
erties of the South Pole ice and estimate the absorption due to solid impurities 
present in the ice. In the discussion we compare our results for South Pole ice with 



measurements made on less pure ice |6L 11, |12|, on extremely pure water |13j, and 



on LiF |n| as an example of a highly transparent solid. 



4.1 Intrinsic absorption in highly pure ice 

Our discussion is based on features shown in Figs. 3 to 5. We focus attention 
on three regions of the approximately V-shaped spectra with their minima in the 
visible region. On the ultraviolet side the absorption in both ice || [ll], [l^] and 



ionic crystals [17] decreases exponentially at wavelengths somewhat longer than that 
corresponding to the electronic band gap energy (~ 8 eV for ice, ~ 13 eV for LiF). 
The slope of this so-called Urbach tail is believed to be governed by exciton-phonon 
interactions; for LiF it becomes more shallow with increasing temperature 17]. 

On the long wavelength side of the minimum, absorption rises approximately 
exponentially as a function of wavelength, but the window of transparency is much 
narrower for ice than for ionic crystals. Ionic crystals are transparent throughout the 
entire visible region, and their exponential absorption spectrum in the infrared is the 
result of multiphonon processes, in which a photon is absorbed with the emission of 
n phonons. Phenomenological models |l8j lead naturally to an exponential depen- 
dence of the absorption coefficient on wavelength, in agreement with measurement 
over six orders of magnitude, for n up to at least seven phonons. The argument 
in the case of ionic crystals is that the contributions of overtones are generated 
by processes that depend on some coupling constant, with the nth overtone going 
as (coupling constant)™ -1 . This leads naturally to an exponential dependence on 
overtone number. 

In contrast, ice is a molecular crystal in which the U-jjO molecules are connected 
by weak OH — O hydrogen bonds. Each O atom in the lattice is coordinated to 
four H atoms, forming strong intramolecular bonds with two of them and hydrogen 
bonds with the other two. Ice absorbs much more strongly than do ionic crystals in 
the red region of the visible. The reason is that ice has additional degrees of free- 
dom associated with the stretching and bending of H-O-H bonds in individual water 
molecules. These internal modes lead to absorption in the red region. The funda- 
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mental stretching mode peaks at ~ 3 /im. Numerous overtones and combinations of 
stretching and bending modes with each other and with the lattice modes lead to 
a roughly exponential slope from a minimum at about 400 nm to a large peak at 3 
/im. The extension of this exponential to values of the absorption coefficient below 
10 -2 m -1 can be seen only if the ice is very pure and if the measurement technique 
is able to discriminate absorption against scattering. 

By application of the same reasoning that applies to multiphonon absorption, we 
suggest that intramolecular vibrations depend on a coupling constant and that the 
overtones decrease as increasing powers of the coupling constant, thus leading to an 
exponential dependence of absorption on wavelength in the red and infrared. Since 
this conjecture assumes no growth in the number of processes with overtone number, 
the actual rate of increase of absorption with wavelength may, for large overtone 
numbers, be less rapid than exponential. (We will make use of this assumption of 
an exponential absorption in Section 4.5 when we estimate the contribution of dust 
to absorption.) 

Impurities, dislocations, vacancies, and their associated electronic defects con- 
tribute to absorption and scattering. Almost no measurement techniques, with the 
exception of our pulsed laser technique, are designed to be able to distinguish scat- 
tering from absorption. These defects tend to fill in the minimum in the spectrum 
to a degree that depends on their concentration and size distribution. 

Figure 3 compares absorption data and the calculated magnitude of Rayleigh 
scattering for ice (solid curves and open symbols) and water (dashed curves and solid 
points) at wavelengths near the absorption minimum. The most striking feature of 
the data is the large discrepancy between values of the absorption near the minimum. 
For water, the absorption decreases in order from the water of Lake Baikal E^ ], 
to the Pacific Ocean [^(J (the DUMAND point), to "pure sea water" (with no 
particulate content) as summarized by Smith and Baker plfl , to the curve labeled 
"purest water", which gives the results of Quickenden and Irvin |13| , who took 
extraordinary measures to eliminate both inorganic and organic impurities, and 
finally to the point labeled SNO. The SNO result is a recent measurement by Noel 
and Mes [23], who went to great lengths to purify the water to be used in the 
Sudbury Neutrino Observatory experiment. 

The Lake Baikal group has recently used a pulsed laser at 475 nm to measure 
scattering separately from absorption. They find values of A sca t ~ 8 — 23 m , 
with {cos 9) ~ 17° [24|. None of the authors of the other studies on water made 
a case for having been able to distinguish Rayleigh scattering from absorption. In 
fact, the data for purest water approach asymptotically (to within a factor two) the 
calculated Rayleigh scattering line for pure water. 

The measurements in Fig. 3 labeled "lab ice" represent the contributions of a 
number of researchers [11, O], none of whom took steps to measure scattering. 
In contrast, the AMANDA experiment measured both the scattered and absorbed 
components separately for ice that is probably the purest naturally-occurring solid 
material on earth. (The concentration of impurities in South Pole ice is discussed 
in Ref. |22|.) At the minimum in the curve for lab ice, the absorption measured by 
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AMANDA is about an order of magnitude lower than the Grenfell and Perovich [|ll|] 
data. Fig. 4 shows the region in more detail, with A -1 as the abscissa on a linear 
scale so that the exponential dependence can be more clearly noted. 

The reader may have noticed that, in the interval ~ 0.2 to ~ 0.45 /im, the data 
for lab ice in Fig. 3 have a slope approximating the Rayleigh A -4 law, but about 
two orders of magnitude higher than the calculated contribution of Rayleigh (or 
thermal) scattering for pure ice indicated by the solid straight line. In Section 4.2 
and in Q we suggest that the apparent absorption of lab ice in this interval may be 
due to Rayleigh scattering from defects in lab ice. The dotted line in Figs. 3 and 4 



is an exponential fit to the data of Grenfell and Perovich [11] in the interval 0.6 fim 
to 0.8 /im. It represents our conjecture as to the intrinsic absorption of ice in the 
red region in the absence of impurities. 

Fig. 5 compares the absorption in ice and a number of transparent solids over 
a broad range of wavelengths. We draw attention to the lines representing fourth- 
power laws, A -4 , which suggest scattering from crystalline defects with dimensions 
much smaller than the wavelength. The techniques used to obtain these data could 
not have distinguished scattering from absorption. Thus, it is tempting to conclude 
that the portions of the absorption curves that exhibit the A -4 dependence are 
mainly due to impurity scattering instead of being due to absorption. In [Q we 
estimate the concentration, size, and type of defects necessary to account for the 
A -4 fits to the data. 

4.2 Absorption by solid dust particles in ice 

We now focus on the AMANDA absorption data [|| taken at depths of 830 m (shown 
in Figs. 3 and 4 as triangular symbols) and at 970 m (shown in Figs. 3 and 4 as square 
symbols). Additional data, taken at a depth of 910 m, were indistinguishable from 
those taken at 830 m. The data || show convincingly that the absorption deviates 
more from the exponential extrapolation at 970 m than at 830 and 910 m. Our 
interpretation is that, in addition to the intrinsic absorption by perfect ice, dust 
particles with a greater concentration at 970 than at 830 and 910 m absorb light 
with a coefficient that decreases with increasing wavelength. This interpretation is 
strongly supported by our model of age vs depth of South Pole ice, by means of 
which we predicted dust concentration to increase to a maximum at a depth slightly 



greater than 1000 m (ref. [14]). 

We now derive the absorption coefficient as a function of wavelength for dust 
in South Pole ice, and we determine the ratio of dust concentration at 830 m to 
that at 970 m. When trying to extract absolute values for the contribution of solid 
impurities to the absorption coefficient, one encounters the problem of estimating 
the uncertainty involved in the exponential extrapolation of the intrinsic absorption 
coefficient of ice shown as the dotted line in Figs. 3 and 4. We developed the following 
way to determine the wavelength dependence of the absorption of solid impurities, 
normalized to some particular wavelength, independent of the unknown intrinsic ice 
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absorption. We write the measured absorption coefficient a exp t(h, A) as 

a e x P t(h, A) = a int (X) + a dust (h, A) = a int (X) + c dust (h)a dust (X), (40) 

which exploits the facts that the intrinsic ice absorption coefficient ctint{X) and the 
absorption cross section of solid impurities a dust (X) are independent of depth, and 
the concentration of impurities, Cdust(h), is independent of wavelength. This assumes 
that the nature of solid impurities is independent of depth and only the concentration 
changes^. 

By using eq. (|40| ) for two different wavelengths and two different depths, we can 
eliminate the unknowns ati n t and Cdust(h) and express the ratio of the absorption 
cross section at two wavelengths in terms of only experimentally measured quantities. 
The result is 

&dust(Xl) _ Oexptjhl, Aj) — aexpt(h-2, Aj) , 

<?dust{X2) a exp t(hi,X2) — a exp t(fi2, A2) 
We will use for the analysis only the data for the six wavelengths that give the 
smallest statistical error. For hi we use data from phototubes at 830 m (around 
AMANDA level 3) and for h 2 we use data at 970 m (around AMANDA level 17). 
We normalize all impurity absorption cross sections to adust (475 nm) , which has the 
most statistics of all the distributions. Figure 6 shows the result of this analysis. 
To test the natural hypothesis that the absorbing impurity is atmospheric dust that 
accumulated in the ice by precipitation onto the surface at the time the ice formed 



from snow [15|, we have compared our results with measurements of the absorption 



spectra of various samples of atmospheric dust at U. S. locations [16]. As can be 
seen, our results are consistent with the hypothesis that we do observe absorption 
by dust. However, it must be noted that the model-independent analysis does not 
have small enough error bars to exclude, e.g., a constant absorption coefficient. 

In Fig. 6 we have also included estimates based on the assumption that the 
intrinsic absorption spectrum of dust-free ice follows the exponential shown by a 
dotted line in Fig. 4. The good fit shown in that Figure was obtained by the 
expression 

a int (X) = 81exp(-A /A) cm" 1 , (42) 

with the value Ao ~ 6.7 /iin being determined not only by a fit to the measured 
absorption in the near IR, but also requiring self-consistency of the Ansatz (|40|). 
That is, the values obtained for dust concentration Cd U st(h) using the experimental 
values a exp t(h, A) and the fitted oti n t(X) should not depend on A; neither should the 
derived adust (A) depend on h. 

Using (|40|), we can derive the following expression for the ratio of the dust 
concentration at the lower level of the detector to the upper level: 

(■dust 

(970 m) _ a expt {%l^Xi) (AO 



('dust 

(830 m) a expt (830, Ai) - a int (Xi) 



(43) 



1 Royer et al. inferred from the optical scattering properties of melted ice core samples from 
Dome C that the particle size distribution and complex refractive index of dust show no depth 
dependence, although it is not clear to us that these conclusions were very strongly supported by 
their data. 



14 



We see in Fig. 7 that all wavelengths give the consistent result ~ 1.4 for this ratio. 

Subtracting now the fitted intrinsic absorption coefficient from the measured 
one at each wavelength, we finally obtain values for the absorption length (i.e. the 
inverse of the absorption coefficient) due to dust at the two AMANDA levels. We find 
that the absorption length due to dust increases with wavelength and decreases with 
depth, with typical values being 200-300 meters, as shown in Fig. 8. This is consistent 
with estimates based on measured dust content and composition in South Pole ice. 
It must be kept in mind, though, that the exponential fit to the intrinsic absorption 
coefficient is only a phenomenological procedure. To the statistical error bars in 
Fig. 8 must therefore be added a systematic uncertainty of unknown magnitude. 

We note that our results at 830 m depth give values for the absorption coefficient 
that represent the purest ice ever measured. Since a residual dust component is 
known to exist also at these depths, the results for a a b s must be interpreted as 
upper bounds for the intrinsic absorption coefficient of pure ice. 

For completeness, we give in Appendix A a brief description of scattering on 
dust grains, expected to be the dominant contribution to the turbidity at larger, 
bubble-free depths. 

4.3 A three-component model of absorption in dusty ice 

To summarize the results of this section, Fig. 9 shows a three-component model which 
we expect to apply to absorption by ice with an extremely low content of dissolved 
impurities. Component 1, which dominates in the ultraviolet, is the exponential 
Urbach tail, given by 

a uv = A uv exp(-0.4818 • A), (44) 

where Ajjy is proportional to the density of insoluble dust in ice. As with ionic 
solids, the slope of the Urbach tail may show a weak dependence on measurement 
temperature. 

Component 2, which dominates in the near-UV and blue, is the contribution of 
insoluble dust, which we assume to have a size distribution similar to that of the 
aerosols studied by Lindberg and co-workers |pf| . Fitting their absorption data to 
a power law in the region 300 < A < 700 nm leads to a dependence proportional 
to A -2 for measurements made both in 1974 and 1994. We assume the same A -2 
dependence for dust in South Pole ice: 

adust = A dust \~ 2 . (45) 

Component 3 is the exponential rise in the red and infrared, as parametrized in 
Eq.g|. 

From section 4.2, we take A^ust at a depth of 970 m to be 1.4 times that at 
830m. The family of curves in Fig. 9 is for values of A^ust increasing incrementally 
by factors of 1.4. 

None of the curves in Fig. 9 fits the data for laboratory ice, nor should they, as we 
noted earlier, the contribution of scattering was not distinguished from absorption in 
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those experiments [||, 11, 12]. Even if scattering had been excluded, the laboratory 
ice had been made from not very pure water and thus probably contained dissolved 
impurities which contributed their own absorption bands to the ice. 



5 Discussion and Conclusions 

The pulsed laser method, applied via optical fiber to deep ice, provides a highly suc- 
cessful method of independently measuring absorption and scattering as a function 
of wavelength, over the interval 400 to 650 nm. 

Previous laboratory studies of absorption of ice at visible and near-visible wave- 
lengths where the scattering length is comparable to or less than the absorption 
length are probably in error due to failure to determine the two parameters sepa- 
rately and probably also to a greater concentration of impurities in ice made from 
laboratory water than in Antarctic ice. Our work has shown that the intrinsic 
absorption of pure ice decreases to values less than a few times 10~ 3 m" 1 at wave- 
lengths between 300 and 400 nm. Our pulsed laser results support the conjecture 
that the absorption increases exponentially with wavelength in the red region. We 
think it would be worthwhile for quantum chemists to seek a theoretical basis for 
the exponential behavior. 

The optical measurements by the AMANDA collaboration, initially intended for 
calibration purposes, have revealed surprising properties of pure ice. The remarkably 
low intrinsic absorption in the visible and near UV region may lead to interesting 
applications for physics and astrophysics. Our simple three-component model for 
absorption by pure but dust-bearing ice (Fig. 9) will be subjected to a critical test 
when the data at depths 1500 to 1900 m taken in 1996 become available. In those 
experiments we will measure absorption and scattering at wavelengths of 337, 350, 
and 380 nm as well as in the region 410 to 650 nm. 
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Appendix A. Simulations of scattering on dust grains 

in Antarctic ice 



Although we have reasons to believe that the optical properties of Antarctic ice 
will be dramatically better at depths where the air bubbles have transformed into 
clathrates, for future use we want to estimate the residual scattering on ice crystal 
boundaries, clathrates, dust and other impurities, as well as Rayleigh scattering in 
the ice itself. 

It is only dust (including soot) of these sources that contributes measurably to 
absorption. We start by briefly discussing the scattering properties of dust, given 
the absorption coefficients derived in the previous Section. 

The total concentration of insoluble impurities in recently deposited South Pole 



ice has been estimated to be around 15 ng/g [25|. The size distribution of dust grains 



has also been measured [26, 27] at the South Pole. At another Antarctic location 
(Dome C), deep samples have been taken and analyzed using a scanning electron 
microscope (giving the number vs radius distribution), a Coulter counter (giving 
the volume distribution), and laser nephelometry (giving the combined scattering 
function of a melted ice sample) fllSfl . In the latter work, the data could be well fitted 
by a log normal distribution with a modal radius of 0.25 fim (for the number to log 
radius size distribution). Using this size distribution and applying Mie scattering 
theory, they could get satisfactory agreement with the measured scattering function, 
except at very large angles, where perhaps effects of the non-spherical shapes of dust 
grains appear. 

We have developed Mie scattering programs to verify the findings of Royer et 
al. for the scattering function. However, when simulations for a large detector like 
AMANDA have to be done it is useful to find approximations that give results much 
faster than CPU-time consuming Mie calculations on an event-by-event basis. 

We have found that a surprisingly good approximation to the results in [|U| 
can be obtained by simply using first-order geometric optics (neglecting internal 
reflections) on dust grains, i.e., by just calculating the deviation angle according 
to Snell's law. This is only a function of where on the beam-facing surface of the 
sphere the light ray enters. We generate a radial variable r at random by setting 
r = tJx (x is a random number distributed between and 1 and the square root 
makes the geometric weight for a given r proportional to r). The scattering angle 
for the Monte Carlo simulation can simply be taken as 

56 = 2(arcsin(r) — arcsin(r/n re ;)), (46) 

where n re i is the relative refraction index (~ 1.18) between dust and ice. 

The expression (^) has the drawback that it underestimates large-angle scatter- 
ing somewhat. In contrast to the multiple scattering case on air bubbles discussed 
before, where only r = (cos 6) is of importance, for reasonable possible spacings 
between PMTs in a working detector only a few scatterings on dust will occur and 
therefore the full scattering function is more important. We have thus chosen a 
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strategy which is based on the Henyey-Greenstein approximation to the scattering 
function: 

— - — T- («) 



d{cos6) (l + r 2 -2rcos( 
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This very accurately reproduces the true scattering function once the value of r = 
(cos 0) is known. To find the latter, we calculated it using Mie theory, with the 
known properties of the dust in South Pole ice as input. For the values of Dome C 
|fl5|| , we find r ~ 0.92. For use in a simulation program, Eq. must be integrated 
and inverted. Setting now r = x (x: random number uniformly distributed between 
and 1), we find that a distribution of the form is obtained by choosing 

2r(l- T + r r)(l+xP 

(l-T + 2rr) 2 v ; 

The 1+cos 2 9 distribution applicable to Rayleigh scattering is similarly generated 
by choosing 

cos6 = 2(k% - -k~3 + ~) -1, (49) 

where 

1 1 1 



- + -r + - v / 5- 16r + 16r 2 . (50) 
4 2 8 

There is a difference between the measured dust size distributions at Dome C and 
at the South Pole. In the latter case, the modal radius is around 0.1 /mi (compared 
to 0.25 at Dome C). This makes the scattering distribution markedly wider as a 
consequence of the approach to the Rayleigh scattering regime, which is isotropic 
on the average. 

It is customary to define the albedo u as the ratio of scattering to extinction 
(i.e., scattering plus absorption) cross sections. The absorption in dust depends on 
the imaginary part of the refractive index, which has been measured to be typically 
between —0.003 and —0.007 for atmospheric dust [28|. From the Mie computations 



(using again a log-normal size distribution) we obtain values of the albedo between 
0.9 and 0.96 both at 500 and 350 nm, when the real part of the refractive index of dust 
varies between 1.53 and 1.56 and the imaginary part between —0.003 and —0.007. 
This is not too different from the values inferred from the Dome C data in ref . . 
The values of r = (cos#) are, however, much smaller: ~ 0.70 for 500 nm and ~ 0.74 
for 350 nm (not very much dependent on the value of the imaginary index). The 
results using a so-called Junge law distribution dn/d(logr) oc r c , with c ~ —1.9 are 
similar, except that the albedo depends more strongly on the imaginary part of the 
refractive index. This seems to be due to the effects of Mie-type resonances whose 
exact location in size parameter depends on both the real and imaginary parts of the 
index. The values of the effective scattering power (defined analogously to Eq. (|i~5D) 
are relatively stable and correspond to u> and r both being around 0.8 - 0.9. 

We note that typical scattering lengths on dust grains in the 800 - 1000 m depth 
range are expected to be of the order of 50 m. Therefore, scattering on air bubbles 
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dominates completely. For greater depths, however, the results of this Appendix 
should be of use. 
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Figure Captions 



1. Method for measuring scattering and absorption in deep ice: pulses at a de- 
sired wavelength are sent down one of 80 optical fibers to a diffusing sphere mounted 
immediately underneath a photomultiplier (PMT) module, and the distribution of 
arrival times is measured at one of the neighbouring PMTs. 

2. Distribution of arrival times as measured and as fitted to the analytical 
formula Eq. 13 for an emitter and a receiver separated by 10 m. The results for 
wavelengths 410 and 610 nm are shown as open and filled circles, respectively. The 
fits shown give the same value for the scattering length (around 10 cm) for the 
two cases, but absorption lengths of around 230 m and 10 m, respectively. 

3. Comparison of absorption spectra in the wavelength interval ~ 200 to ~ 800 
nm for ice (solid curves and open symbols) and water (dashed curves and solid 
points). 

4. Absorption as a function of wave number in /urn -1 for laboratory-grown ice 
crystals (open circles) and for South Pole ice at level 3 (830 m; triangular symbols) 
and level 17 (level 970 m; square symbols). The dotted curve is an exponential fit 
to data for lab ice at 600 to 800 nm. 

5. Comparison of absorption spectrum in ice and a number of other transparent 
solids. The straight lines representing an inverse fourth-power law, A -4 , are meant 
to show that portions of each spectrum are due to Rayleigh scattering. 

6. Wavelength dependence of the absorption on dust inferred from the AMANDA 
laser measurements. The normalization is arbitrarily set to unity at A= 475 nm. 
Open squares, and the error bars, refer to the model independent analysis where 
the intrinsic absorption in ice has been subtracted. Triangles and circles are values 
from AMANDA depths 830 and 970 m, respectively, using an exponential fit to 
the intrinsic ice absorption. The band within the dashed lines corresponds to the 
wavelength dependence measured for atmospheric dust at a US location |16[ |. 

7. Ratio of dust concentration at 970 m and 830 m as obtained by the analysis 
in this paper. Note that the ratio does not depend on wavelength, which is a 
consistency check of the scheme employed. 

8. The absorption length (i.e. the inverse of the absorption coefficient) given by 
scattering on impurities in the ice, as a function of wavelength for AMANDA 830 
m level (filled triangles) and 970 m (filled squares). The values were obtained by 
subtracting the best fit exponential for the intrinsic absorption from the measured 
absorption. 

9. Three-component model of absorption by ice containing insoluble dust but no 
dissolved impurities. The concentrations of dust increase by factors of 1.4 for each 
curve. 



23 



Wavelength (nm) 


\83Um f rr ,\ 
*abs l m ) 


\970m f rr ,\ 


a i n 

41U 


91Q Q + 19 ^ 




415 


225.6 ±6.2 


168.7 ±23.6 


420 


227.4 ±4.0 


163.9 ±7.1 


435 


207.4 ± 1.5 


166.3 ± 1.9 


450 


171.3 ± 1.3 


146.0 ±2.1 


475 


117.1 ±0.7 


104.6 ± 1.2 


500 


68.9 ± 1.3 


64.1 ± 1.0 


530 


38.3 ± 1.0 


39.5 ±0.7 


590 


12.0 ±0.1 


10.5 ±0.3 


610 


10.3 ±0.3 





Table 1. Absorption lengths X a b s (i.e. the inverse of the absorption coefficients a a f, s 
defined in Eq. ( |6|) for several laser wavelengths at two different depths, around 
830 and 970 m, determined by the AMANDA calibration system in the South 
Pole ice. For 410 and 610 nm, only data for 830 m could be obtained. 
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